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Abstract. We analyse the Chandra dataset of the galaxy cluster MS 1008. 1 — 1224 to recover an estimate of the gravitating 
mass as function of the radius and compare these results with the weak lensing reconstruction of the mass distribution obtained 
from deep FORS 1 -VLT multicolor imaging. Even though the X-ray morphology is disturbed with a significant excess in the 
northern direction suggesting that the cluster is not in a relaxed state, we are able to match the two mass profiles both in absolute 
value and in shape within la uncertainty and up to 1 100 hg kpc. The recovered X-ray mass estimate does not change by using 
either the azimuthally averaged gas density and temperature profiles or the results obtained in the northern sector alone where 
the signal-to-noise ratio is higher. 
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1. Introduction 

As the largest virialized objects in the Universe, galaxy clusters 
are a powerful cosmological tool once their mass distribution is 
univocally determined. In the recent past, there have been sev- 
eral claims that cluster masses obtained from X-ray analyses of 
the intracluster plasma, taken to be in hydrostatic equilibrium 
with the gravitational potential well, are significantly smaller 
(up to a factor of two; but see Wu et al. 1998; Allen 1998; 
Bohringer et al. 1998; Allen et al. 2001) than the ones derived 
from gravitational lensing (see Mellier 1999 for a review). 

In this paper we report on the mass distribution of the clus- 
ter MS 1008. 1 — 1224 by combining the results from weak lens- 
ing analysis of deep FORS 1 - VLT images with those obtained 
from a spatially-resolved spectroscopic X-ray analysis of a 
Chandra observation. MS 1008. 1 — 1224 is a rich galaxy clus- 
ter at redshift 0.302 that has been part of the Einstein Medium 
Sensitivity Survey sample (Gioia & Luppino 1994) and of the 
CNOC survey (Carlberg et al. 1996). Lombardi et al. (2000) 
presented a detailed weak lensing analysis of the FORS 1 -VLT 
data. Figure [l] shows the X-ray isophotes overplotted to the op- 
tical V-band image. In the following we adopt the conversion 
1 arcmin = 376 kpc (z = 0.302, H = 50 h 5Q km s _1 Mpc" 1 , 
Cl m = 1 - Q a 
fidence level). 




0.3) and quote all the errors at la (68.3% con- 



35= 30° 
Right Ascension (2000) 



Fig. 1. An adaptive smoothed exposure-corrected image in the 
Chandra [0.5-2] keV band is overplotted to the R-band image 
(the white areas in the upper end show masked regions). The 
contours are spaced with steps of log 2 from the minimum of 
2.1 x 10~ 9 photon sec" 1 cm -2 arcsec -2 . 
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Fig. 2. Residuals above 2 in the X-ray brightness distribu- 
tion. Given two regions, A and B, of (length, width) = (8 
arcsec, 20°) located symmetrically with respect to the center 
and with Ca and Cb observed counts respectively, the resid- 
uals are estimated as a = (Ca — Cb)/\/Ca + Cb- An az- 
imuthal scan with step of 20°/3 was done to smooth the map. 
The most significant excess is along 63° ± 10° (anticlockwise, 
from X-axis) with a power of 0.87 estimated as the fraction of 
regions with higher excess along a fixed direction. 



2. X-ray mass 

We retrieved the primary and secondary data products from 
the Chandra archive. The exposure of MS 1008. 1 — 1224 was 
done on June 11, 2000 using the ACIS-I configuration. We 
reprocessed the level=l events file in the Very Faint Mode 
and, then, with the CtiCorrector software (v. 1.38; Townsley 
et al. 2000). The light curve was checked for high background 
flares that were not detected. About 44.0 ksec (out of 44.2 ksec, 
the nominal exposure time) were used and a total number of 
counts of about 20 000 were collected from the region of in- 
terest in the 0.5-7 keV band. We used CIAO (v. 2.2; Elvis 
et al. 2002, in prep.) and our own IDL routines to prepare 
the data to the imaging and spectral studies. The X-ray cen- 
ter was fixed to the peak of the projected mass from weak 
lensing analysis (Lombardi et al. 2000) at (RA, Dec; 2000) = 
(10 h 10 m 32?68, -12°39'58.8"). Note that the maximum value 
in a 5"-smoothed image of the cluster X-ray emission is at (RA, 
Dec) = (10 h 10 m 32!44,-12°39'55.6"), i.e. less than 5 arc- 
sec apart from the adopted center. With respect to the adopted 
center, a clear asymmetry in the surface brightness distribution 
is however detected, suggesting an excess in emission in the 
northern region (see Fig. 0). 



Fig. 3. Map of the distribution of the gas temperature accord- 
ing to the structures present in the X-ray surface brightness (see 
Fig. [|and Sect.2.1). The code color represents also the varia- 
tion in the range (—3.1, +2.1)ct with respect to the global value 
of 7.2 keV. 



2. 1. X-ray analysis 

We detected extended emission at 2a confidence level up to 4. 1 
arcmin (1.55 Mpc) and we were able to extract a total of four 
spectra with about 2 000 source counts each to put reasonable 
constraints on the plasma temperature profile up to 1140 kpc. 
The contribution from the source to the total count rate de- 
creases from more 95% in the innermost spectrum to about 
30% in the outermost one. We also evaluated the distribution 
of the plasma temperature values in sectors according to the 
asymmetrical surface brightness shown in Fig. ||. These best- 
fit values for a two-dimensional map and for an azimuthally 
averaged profile are represented in Figs. || and|j. 

We obtained the Redistribution Matrix Files (RMFs) and 
Auxiliary Response Files (ARFs) by using the CIAO rou- 
tines mkrmf and mkarf with the QEU files included in the 
CtiCorrector package. An emission from an optically- 
thin plasma (Mekal - Kaastra 1992, Liedhal et al. 1995, in 
XSPEC v. 11.1.0 - Arnaud 1996) with the metal abundance 
fixed to 0.3 times the solar value (Anders & Grevesse 1989) 
and absorbed from the interstellar medium parametrized us- 
ing the Tubingen-Boulder model (tbabs in XSPEC; Wilms, 
Allen & McCray 2000) was adopted to reproduce the observed 
spectra. A galactic column density fixed to 7.0 x 10 20 cm -2 
(from radio HI maps in Dickey & Lockman 1990) was as- 
sumed. A local background was adopted also considering the 
relatively high column density of this field with respect to the 
blank field available for the same CCD and the proper obser- 
vational period. The overall spectral fit of the counts collected 
within 1100 kpc from the adopted center provides an emission 
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Fig. 4. Projected (triangles, dotted line) and deprojected (diamonds) gas temperature profiles. The solid line shows the best-fit 
applying the hydrostatic equation in combination with the deprojected electron density as discussed in Sect. 2.2. The dashed line 
indicates the best-fit from a polytropic model. Results for the North excess (left) and for azimuthally averaged profile (right) are 
shown with overplotted the ler range of emission- weighted temperature up to 110 kpc. 



weighted temperature of 7.2+J' g keV and a bolometric lumi- 
nosity of 1.6 x 10 45 erg s" 1 (0.9 x 10 45 in the 2-10 keV band). 



2.2. X-ray mass profile 

In accordance with the weak lensing analysis described in the 
following section, we assume a spherical geometry for both the 
dark matter halo and the X-ray emitting plasma (note that neg- 
ligible effects, when compared with our statistical uncertain- 
ties, can be introduced on the mass estimate due to the aspher- 
ical X-ray emission, see, e.g., Piffaretti et al. 2002). The values 
of gas density and temperature in volume shells are recovered 
from the projected spectral results as described in Ettori et al. 
(2002). To measure the total gravitating mass M to t, we then 
constrained the parameters of an assumed mass model by fit- 
ting the deprojected gas temperature (shown in Fig. with the 
temperature profile obtained by inversion of the equation of the 
hydrostatic equilibrium between the dark matter potential and 
the intracluster plasma, 



G/im F 



3 M to t,model(< r) _ d(n c kT) 



dr 



(1) 



In this equation, /i — 0.6 is the mean molecular weight in 
a.m.u., G is the gravitational constant, m p is the proton mass, 
and n c is the deprojected electron density. We considered both 
the King approximation to the isothermal sphere (King 1962) 
and a Navarro, Frenk & White (1997) dark matter density pro- 
file as mass models (see details in Ettori et al. 2002). By fit- 
ting the temperature profile in Fig. || we measured the best fit 
parameters (r s ,c) = (646 ± 390,4.2 ± 0.9) for a King and 
(1122 ± 287, 2.3 ± 0.8) for a NFW mass model. 



Using a /3-model (Cavaliere & Fusco-Femiano 1976) and 
the ROSAT HRI surface brightness profile detected up to about 
765 kpc at the 2er level, Lewis et al. (1999) measured A/ to t = 
(3.8±0.6) x 10 14 M Q , assuming an isothermal gas temperature 
of 7.3 keV. If we apply a /3-model to our surface brightness 
profile and a polytropic function to the gas temperature pro- 
file, we obtain r c = 0.264±° °:| Mpc, (3 = 0.618±g;g§f and 
7 = 0.96^q']J. The derived mass estimate is lower by 10-20% 
(by 12% at 765 kpc) than the one in Lewis et al. (1999). The 
best-fit mass models give M to t,King = (4.5 ± 1.2) x 10 14 M Q 
and M to t, nfw = (3.2 ± 0.5) x 10 14 M Q , in agreement with 
the results obtained from each independent sector of the tem- 
perature map in Fig. |] [e.g., the region to North, which has 
a higher signal-to-noise ratio due to the excess in brightness, 
gives Mtot.Kmg = (4.5 ± 0.8) x 10 14 M Q and M tot , NFW = 
(4.1 ±0.6) x 1O 14 M ]. 



3. Weak lensing mass 

A weak lensing analysis of MS 1008. 1-1224 using VLT- 
FORS 1 images was carried out by Lombardi et al. (2000) and 
is summarized below. A parallel weak lensing analysis carried 
out by Athreya et al. (2002), leads to a mass estimates in agree- 
ment with the one presented here. 



3. 1. Weak lensing analysis 

The weak lensing analysis was performed independently on 
the four (B, V, R, and I) FORS-1 optical images using the 
IMCAT package (Kaiser et al. 1995; see also Kaiser & Squires 
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Fig. 5. Projected X-ray total gravitating mass profile (with relative error) obtained directly from the electron density and tem- 
perature profiles is shown as histogram. Cumulative projected mass profiles from X-ray best-fit mass model (King: dashed line; 
NFW: dot-dashed line) and weak lensing analysis (squares/B band, diamondsN band) are overplotted. The deviations in a and 
in the value of the ratio between X-ray best-fit mass model and weak lensing projected masses are plotted in the panels at the 
bottom. 



1993). After generating source catalogs, we separated stars 
from galaxies; the measured sizes and ellipticities of stars were 
used to correct the galaxy ellipticities for the PSF (see Kaiser 
et al. 1995; Luppino and Kaiser 1997). The observed galaxies 
were classified as foreground, background, and cluster mem- 
bers using, when available, the redshift from the CNOC survey 
(Yee et al. 1998), or the observed colors otherwise. Finally, we 
used the fiducial background galaxies to obtain the shear field 
of the cluster and, from this, the projected mass distribution 
(see Lombardi & Bertin 1999). 

Our study differs on a few points with respect to similar 
weak lensing analyses: 



- We decided to use a robust, median estimator to obtain 
the local shear map from the galaxy ellipticities instead of 
the more common simple average. In particular, we esti- 
mated the local shear on each point of the map by taking 
a (weighted) median on the observed ellipticities of angu- 
larly close background galaxies (see Sect. 4.4 of Lombardi 
et al. 2000). This way we make sure that a few galaxies 
with poorly determined ellipticities do not significantly af- 
fect the shear estimation. 

- We estimated the background sources redshift distribu- 
tion by resampling the catalog of photometric redshifts of 
Fernandez-Soto et al. (1999) in the Hubble Deep Field. 

- We took advantage of the multi-band observations by per- 
forming the weak lensing analyses on each band separately. 



3.2. Weak lensing mass profile 

The reconstructed two-dimensional mass distribution appears 
to be centered on the cD galaxy and shows elliptical profiles 
oriented in direction North-South. No substructure on scales 
larger than 30" was detected. In order to remove the mass- 
sheet degeneracy (see, e.g., Kaiser & Squires 1993), we fit- 
ted the mass profiles with non-singular isothermal sphere mod- 
els; from this we obtain, for example, M(r < 1 h§Q Mpc) = 
5.3 x 10 14 /j^Mq. We note that, because of the smoothing 
operated in the two-dimensional mass maps, the weak lensing 
mass profile for r < 1' is bound to be underestimate. On the 
other hand, at large radii (say, r > 3'), an error on the removal 
of the mass-sheet degeneracy can in principle lead to an unre- 
liable "total" weak lensing mass estimate. Note that the four 
profiles (from B, V, R and I optical images) agree very well to 
each other, which strongly support the results of our analysis. 

4. Conclusions 

The differential X-ray best-fit mass model has been weighted 
by the relative portion of the shell observed in each ring and, 
then, cumulated up to the radius of 1020 kpc (the outer ra- 
dius of our spatially resolved spectroscopy). In Fig. ^, we plot 
and compare the projected mass profiles of the galaxy clus- 
ter MS 1008. 1-1224 obtained from both the spatially resolved 
spectral analysis of the Chandra observation and the weak lens- 
ing analysis of FORS1-VLT multicolor imaging. The two in- 
dependently reconstructed mass profiles agree very well within 
lcr uncertainty both in absolute values and in the overall shape 
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of the profile. Note that the mass center is fixed to the peak of 
the lensing map density that is consistent with the X-ray peak as 
discussed in Sect. 2. This result does not change when we con- 
sider the different density and temperature profiles observed in 
the northern region where a significant surface brightness ex- 
cess is located and a higher signal-to-noise ratio is available, 
arguing for the robustness of the X-ray mass estimates once gas 
density and temperature distributions can be properly mapped. 
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